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The primary objectives of Experiment 74-42/2, which flew as part of the 
SPAR VIII payload on November 18, 1980, were melting and resolidification of 
a glass specimen in a containerless condition and the retrieval and examination 
of the specimen. The absence of container contact was to have been assured 
by use of a single-axis acoustic levitation system. However, the sample con- 
tacted a wire cage after being held without container contact by the acoustic 
field for only approximately 87 seconds. At this time, the sample was still 
molten and, therefore, flowed around the wire and continued to adhere to it. 

This report concentrates on an analysis of why the sample did not remain levi- 
tated free of container contact. The experiment is described, experimental 
observations are discussed and analyzed, and recommendations for future single- 
axis acoustic levitation experimentation are provided. 
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TECHNICAL MEMORANDUM 


ANALYSIS OF SPAR VI 1 1 SINGLE-AXIS 
LEVITATION EXPERIMENT 

I. EXPERIMENT DESCRIPTION 


A. Objectives 

Experiment 74-42/2, which flew as a part of the Space Processing Applications 
Rocket (SPAR) VIII payload on November 18, 1980, had as primary objectives the 
melting and re solidification of a glass specimen in a containerless condition, and 
the retrieval and examination of the specimen. The absence of container contact 
was to have been assured by use of a single-axis .-caustic levitation system designed 
and built by Intersonics, Inc. The experimental objedtives were degraded by two 
events. First, the sample contacted a wire cage after being held without container 
contact by the acoustic field for only approximately $7 sec. At this time, the 
sample was still molten and therefore flowed around the wire and continued to 
adhere to it. Second, the sample was lost (and probably destroyed) because of 
failure of a parachute and the resulting crash of the SPAR payload. We will 
concentrate here on analysis aimed at the first problem; that is, why the sample 
did not remain levitated free of container contact. 


B . Apparatus 

The single-axis acoustic levitation system consists essentially of an acoustic 
generator and amplifier which drive linear oscillations of a cylinder. The circular 
area forming the cylinder end oscillates in a direction parallel to the cylinder axis, 
producing sound waves which are reflected by a circular plate which is co-axial 
with the driving cylinder (Fig. 1) . The reflected sound wave interferes with the 
primary wave, producing a three-dimensional pressure field which ideally has 
cylindrical symmetry about the cylinder axis. Thi^ pressure field has minima at 
nodal points along the cylinder axis which correspond to the bottoms of three- 
dimensional potential wells. An object placed in the field will move to the center 
of the well or oscillate around or through it depending on the initial momentum and 
on the absence of other dominant forces. Figure 2 shows qualitatively the shape 
of these potential wells. The potential is generally considered very nearly of the 

form l/2Kr for the radial dependence and l/2Kz for the axial dependence. That 
is, in this approximation, the motion of a levitated object (neglecting other forces) 
can be described in terms of superposition of two independent motions. One of 
these is a simple harmonic oscillation parallel to the z-direction. 9 he other is motion 
due to a central force of magnitude Kr. This motion can be straight line, circular, 
or elliptical depending on initial conditions. The details of these motions and the 
nature of deviations from them will be discussed in depth in this report. 

The configuration for Experiment 74-42/2 is shown in Figure 3. In addition 
to the sound source and reflector which produce the levitating field, there are 
heater rods around the volume occupied by the sample which provide the heat 
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Figure 1. Acoustic levitator configuration showing extent of sample motion. 
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Figure 2. Acoustic potential well (radial, with shape factor b) 
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Figure 3. Acoustic levitator configuration for SPAR VIII. 

necessary to nielt the sample. There is an injection system consisting of a springloaded 
cage designed to hold the sample motionless against the reflector for protection during 
launch. At the beginning of the experiment, the cage is moved along the z-axis so 
the sample is free to move without contact with any container. After processing, the 
cage is designed to retract (moving in the +z-direction) so that the sample is recap- 
tured and held motionless during re-entry and impact. There is an optical system, 
consisting of a mirror and a motion picture camera, for observation of the sample during 
the course of the experiment. 

The experiment package is mounted on the front end of the SPAR payload (Fig. 4). 
The entire 74-42/2 package is within less than 2 of the payload center of mass. (The 
payload center of mass is designed to be on the z-axis and near the geometric center 
of the payload.) The orientat: in of the experiment is shown in Figure 3. The axis of 
the levitator (geometric axis through the center of driver and reflector) is parallel to 
the spacecraft z-axis. The camera mounting arrangement produces a view in which the 
positive y-direction is up and the positive ^-direction is to the right. That is, 
although the image produced of the sample is of the side toward the positive x- direc- 
tion , the motion observed is that which would be seen by looking along the x-axis 
from - to + with +y upward. In addition, Intersonics, Inc., reported that the camera 
was rotated by approximately 6 degrees in its mounting so that in image, the levitator 
acoustic axis would appear to be rotated counterclockwise by approximately 6 degrees. 
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Figure 4. SPAR R-19 science payload May 19, 1980, configuration. 


C . Sample 

The sample processed in the Intersonics levitator and furnace system was 
reported by Dr. Ralph Happe of Rockwell (the Principal Investigator) to be a 
gallia(41%)-calcia(36%)-silica(23%) sphere having a mass of 0.595 gm. Measurements 
of neither diameter nor density were available for the flight sample. However, using 
the molar fraction concentrations given by Dr. Happe and the Handbook of Chemistry 
and Physics (50th Edition) values for densities of gallia, calcia, and silica, a density 

of 4.43 gm/cm was computed (see Appendix A). This yields a value of the sphere 
diameter of 0.635 cm (0.250 in.). 


D. Experiment Timeline 

Figure 5 shows the planned experiment timeline superimposed on the SPAR 
trajectory. The experiment functioned as planned until the sample contacted the 
wire cage after approximately 87 sec of experiment time (launch +204 sec). 


E. Data Collection 

The SPAR acoustic single-axis levitation experiments yielded data of four 
separate types. First, a 16 mm motion picture camera recorded visual observations 
of the motion of the sample sphere in the levitator from sphere injection throughout 
the heating and cooling phases of the experiment. Second, a three-axis acceler- 
ometer (on SPAR VIII) allowed monitoring of linear accelerations at the accelerometer 
positions throughout the low-g phase of the flight. Third, a rate gyro system 
provided continuous monitoring of angular rates about roll, pitch, and yaw axes. 
Fourth, housekeeping data for the experiment were provided continuously. These 
included measurements of acoustic driver voltage, current, and phase as well as 
vibr .tor amplitude and temperature. 
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Figure 5. SPAR project NIKE -BLACK Brant VC rocket 
mission profile, events and conditions. 


II. EXPERIMENTAL OBSERVATIONS 


A. Primary Data 


1. Film Record. 

The motion picture record of the experiment shows the sphere in the cage 
assembly prior to injection into the levitating acoustic field. Injection is observed 
at the expected time. The sphere then is seen to develop a rotational motion about 
an axis internal to the sphere. This rotation appears to build up to a rate which 
becomes comparable to the camera frame rate (24 frames per sec). This fact, 
coupled with low contrast between sample and background , makes the observation of 
sample spin rather qualitative. 

The sample is also seen to begin an oscillatory motion which has components 
in the direction of the levitator axis and along a levitator radius. Since only one 
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camera view is available, we are able to observe only a projection of the sample 
motion onto the camera viewing plane. This oscillatory motion is observed to persist 
for approximately 87 sec after sample injection. Drifts are observed in the apparent 
centers of motion both axially and radially during this 87-sec period. A; the end 
of this period the sample is observed to move abruptly a short distance (about a 
sample radius) and then come to rest, apparently adhering to the sample cage. The 
sample was molten at that time and appeared to have flowed around one of the cage 
wires and centered itself on the wire. 

2. Accelerometer Data. 


The three -axis accelerometer records are shown in Figures 6 through 8 . The 
A02 record indicates an oscillatory behavior of that component of acceleration having 

- 5 

an amplitude no greater than approximately 1.5 x 10 g. These oscillations are 

' 6 

about a steady background level of approximately -5 x 10 g. Assuming that the 
zero calibration of the uccelerometer represents zero gravity, the excursions of 

- 5 

measured acceleration levels for A02 lie between approximately 1 x 10 g and -2 x 

- 5 - 5 

10 g. It is likely that the measured steady component at -5 x 10 g represents 

a zero reading from the instrument. This would imply that, acceleration levels 

- 5 - 5 

fluctuated between approximately +1.5 x 10 g and -1.5 x 10 g. Whichever inter- 
pretation represents the real situation , the magnitude of observed acceleration levels 

- 5 

for A02 never exceeds approximately 2 x 10 g. Accelerometer A03 recorded a 

- 5 

mean background acceleration level of approximately -1.5 x 10 g wit: maximum 

~ 5 

fluctuations about this level of approximately i2.5 x 10 g. Assuming agr ; n that 
the accelerometer zero level is unshifted, the range of acceleration excursions in 

C c _ c 

this direction is approximately +1.0 x 10 g to -4.0 x 10 g. If the -1.5 x 10 g 

reading indicates the ze-o acceleration level, then the range is +2.5 x 10 ^ g to 
-5 

-2.5 x 10 g. Accelerometer A04 indicates a background acceleration level which 

-5 -5 

decays from approximately 1.0 x 10 g to approximately - 0.2 x 10 g over the 
low-g phase of the flight. Maximum excursions from this mean are approximately 
-5 

i - 5 x 10 g. Assuming zero shift in accelerometer calibration, accelerations in this 

-5 -5 

direction lie within the limits +3.0 x 10 g to -2.5 x 10 g. 

3. Angular Rate Data. 

Data are available for SPAR VIII from the payload's rate gyro system giving 
rotation rates about roll, pitch, and yaw axes (Figs. 9 through 11). Rotations 
about each of these yaw axes can potentially give rise to accelerations of the sample 
in at least two different ways. First, centrifugal accelerations can arise due to the 
rate of change of angular position according to 

v fc r = a . 

c 

This can be broken down into motion about each of these independent axes (yaw, 
pitch, and roll). We will have then: 
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Figure 9. Roll rate. 
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where ui , and are rotations about the pitch, roil, and yaw axes respectively; 

and r are distances from the pitch, roll, and yaw axes to the experiment 
ye 
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pe re 
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pe re 
at the experiment location, 
tion axes. 


.. are acceleration components experienced by a sample 
These accelerations are directed outward from the rota- 


In addition , there are 
rates. These accelerations 


accelerations which arise due to changes in the rotation 
are of the form : 


b = ui x r. 


Components of this acceleration will have the following magnitudes: 
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where Cn . i'o , *.i are the time rates of change of rotation rate about the roll, 
pry H 


pitch, and yaw axes respectively (Figs. 12 through 14): and b , t> re » b are 

the components of acceleration experienced by a sample at the experiment location 
due to rates of change of rotation about the respective axes. The accelerations here 
are directed perpendicularly to the rotation axis and to the radius from the given 
axis out to the experiment. 


The flight data must be examined to determine acceleration levels induced by 
residual spacecraft motion. The levitated sample was lost after less than 100 sec of 
experiment time, which fell within the period t = 104 sec to t = 204 sec on the 
experiment timeline. Good angular data exist for this period, so the required 
analysis is possible. During this period, roll rates began at 0.14°/sec, went 
through 0, and climbed to +0.05°/sec. The maximum value (absolute value) of roll 
rate was 0.14°/sec. The average rate of change of roll rate over this period was 
2 

0.0021°/sec . Pitch rates ranged from a beginning value of -0.07°/sec to 0°/sec at 
204 sec. Maximum absolute value was 0.07°/sec. Average rate of change of pitch 
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Figure 14. Rate of change of yaw rate. 

2 

rate was 0.0007°/sec . The yaw rate ranged from an initial -0.07°/sec to a final 
rate (at t = 204 sec) of -0.2°/sec. Maximum absolute value of yaw rate was 0.2°/sec. 

2 

Average rate of change of yaw rate was 0.0013°/sec . None of the angular rate 
changes varied greatly from linearity with time. 

Experiment distances from the roll axis are on the order of 0.1 m. Distances 
from the yaw and pitch axes are less than 2 m. These values can be used to place 
upper limits on the magnitudes of accelerations experienced by a sample at the 

experiment location due to and u>. 

a ^ 5.8 * io ' 8 g 

re & 

a = 2.9 x 10 ^ g 

pe 6 

a = 2.5 x io " 6 g. 

ye 6 


a and a are in the same direction . but their combined effect is still on the order 
re _ ye 

of 10'° g. 
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These acceleration levels are all too small by at least two orders of magnitude 
to have any significant effect on the motion of the levitated sample. 

4. Levitator data. 

Figures 15 through 19 show levitator housekeeping data. These data indicate 
that the levitation system acoustic output was normal throughout the experiment and 
that the vibration temperature was within allowable limits. It can be inferred from 
this that electrically and mechanically the acoustic driver and vibrator functioned 
normally. The loss of sample by drifting into the cage then was caused by some 
other factor or factors. 
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Figure 17. Acoustic driver phase. 
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Figure 18. Vibrator amplitude. 
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Figure 19. Vibrator temperature. 
B . Secondary Data 


1. Digitized position data . 

The motion picture film of the SPAR VIII experiment was analyzed frame by 
frame to obtain a record of the projection of the sample trajectory onto the film 
plane. This set of measurements from the film was done by Intersonics, Inc. The 
data were compiled in tabular form showing frame number and sample coordinates 
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(right, left, top, and bottom sumple edge positions were shown). This tabulated 
data set was provided to the Marshall Space Flight Cenler Space Sciences Laboratory 
for analysis. Here the data were transferred to magnetic tape so that data presen- 
tation and analysis could be more readily accomplished. All the position data dis- 
cussed later were plotted from this tabulated data set. 

In the simplest approximation to the sample motion, the two-dimensional pro- 
jection of the three-dimensional sample trajectory provides a fairly complete descrip- 
tion of the sample motion. As noted previously, the motion in this approximation 
can be described by an axial (z-direction) simple harmonic motion superposed on a 

2 

central force motion [with V(r) = 1/2 Kr ]. The radial motion can further be 
described by a superposition of two linear harmonic oscillators with the same force 
constant K. A measurement of the period of any projection of the radial motion onto 
a line is then sufficient to calculate the force constant. Details of the actual 
trajectory, i.e. , x(t) and y(t), are not interesting in this approximation for this 
application and would require knowledge of x(t ) and y(t ), where t is any value of 
time (initial condition). ° c 

It was observed in the analysis of SPAR VIII experimental results that the 
motion of the sample in the acoustic field does not fo. >w exactly this simple harmonic 
oscillator description. The deviations from t ids model are both interesting and 
critical to the outcome of the experiment. The use of a single- view camera system 
is adequate when the motion is near simple harmonic, but information about sample 
behavior (hence about the nature of the acoustic field) is lost when deviations from 
this motion are significant. The use of an optical system which provides projections 
of the motion on two perpendicular planes would permit a frame by frame plot of 
sample trajectory which would allow calculations of acoustic forces (magnitude and 
direction) as a function of time. 

2. Axial motion . 

Figure 20 shows a plot of the projection of the sample motion onto the z-axis. 
The motion is seen to be cyclic with a varying amplitude and period and with a net 
drift in the mean sample position. Figure 21 shows the amplitude of the axial oscil- 
lation as a function of time. Time variation of the period of the axial oscillation over 
the duration of the experiment is shown in Figure 22. A plot of average axial posi- 
tion with time gives a direct indication of drift of the axial "center of motion" or 
"equilibrium position". This is shown in Figure 23. 

Using simple harmonic oscillator approximation to the axial motion, an approx- 
imate force consta nt c an be extracted from the sample period as a function of time 
(from <o= 2ir/T = vdt/m) . Axial force constant, in this approximation, is shown in 
Figure 24. As indicated elsewhere, this approximation is not adequate for a full 
description of sample behavior, but does serve to provide a lowest order description. 

3. Radial motion . 

Sample radial motion -that is, the projection of motion in the x-y plane onto a 
line in that plane (perpendicular to the camera viewing direction) -is shown in 
Figure 25. The radial motion shows some of the characteristics of the axial motion. 
That is, there is an oscillatory motion which varies in amplitude and frequency over 
the duration of the experiment; and there is a drift in the average sample position. 
The amplitude of the radial oscillations as a function of time is shown in Figure 26. 
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Figure 20. Axial motion of sample center. 
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Figure 21. Amplitude of axial oscillations as a function of time. 
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Figure 25. Radial motion of sample center. 
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Figure 27 shows the time behavior of the period of radial oscillation. "Center of 
motion" drift, or average radial position as a function of time, is shown in Figure 28. 

2 

The 1/Zk r r potential approximation of the radial motion allows calculation of a 

force constant from the oscillation frequency. It is computed as in the i-xial case 
and is shown plotted as a function of time in Figure 29. 

4. Sample trajectory in the y-z plane . 

Appendix B consists of a sequence of plots of sample trajectory projected c~ to 
the y-z plane (camera film plane). Plots are shown in their correct vime sequence, 
with each plot representing luO film frames or approximately 4.17 sec. 


III. OBSERVED PHENOMENA REQUIRING EXPLANATION 


The data presented previously comprise the complete :et of raw datr from the 
flight experiment and a set ot secondary data based ,n sample position as a function 
of time. An aim of this work is to explain the phenomeha observed in the data as 
completely as possible. Key phenomena are the following: 

1) Oscillatory motion of the sample 

a) axial component of oscillation 

b) radial component of oscillation 

2) Rotation of sample 

a) Buildup of rotation abcut internal axis 

b) orientation of rotational axis 

c) magnitude of observed rotation r.ste 

3) Variations in oscillation frequency 

a) axial variations 

b) radial variations 

4) Variations in oscillation amplitude 
a) axial variations 

b> radial variations 
5; Drift in average sample position 

a) axial drift 

b) radial drift 
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6) Sudden motion of sample followed by the sample coming to rest and remain- 
ing at rest for the duration of the experiment. 


IV. FORCES ACTING ON SAMPLE 


Before analyzing sample motion in detail, it is useful to catalogue the forces 
which can act on the sample. These are: 

1) Forces due to the payload and its motion were measured by the low-g 

accelerometers and are derived from the rate gyro data. Measurements indicated 

. 5 

that none of these is larger than a few times 10 g. 

a) Drag from the interaction of the payload with the very tenuous atmo- 
sphere encountered during the low-g trajectory phase of the flight would be expected 

7 

to be less than 10 g. 

b) Impulses due to other experiments or from the payload assembly (e.g.. 
gas venting) were not expected to be a problem, and the accelerometer data yield 
upper limits for these effects. 

c) Centrifugal accelerations at the experiment location due to payload 

6 

rotations were shown to be a few times 10 g or less. 

d) Poincare forces (>o \ r) due to angular accelerations were shown to be 

6 

also of the level of a few times 10 g or less. 

2) The sound field produced by the acoustic driver and reflector has a dc 
component which produces the positioning (or levitating) force. This will be dis- 
cussed in further detail in Section V. This dc field is affected by several factors. 

a) The density of the medium is related to the magnitude of the levita- 
ting force. This implies a temperature dependence. That is, as air temperature 
goes up (and density goes down), the levitating force decreases. 

b) Geometrical effects in the acoustic system can be very important. For 
example, unwanted reflections can produce distortions i.i the axial symmetry of the 
dc sound field. These reflections can be produced by objects such as the sample 
injection assembly which intrude into the sound field. 

c) Changes in moisture content of the levitating medium can result in 
sound absorption wluch can reduce pressure levels, 

3) Acoustic streaming is also discussed in Section V. This flow produces 
aerodynamic forces on the sample. These forces may be important when the sample 
moves very near the cage assembly. 

4) Aerodynamic forces due to sample motion (drag und magnus effect) are 
discussed in Section V. These forces are too small to influence sample motion appre 
ciably . 
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5) Surface tension force h iieeome dominant when the molten sample cent acta 
the cage assembly. This causes the final rapid motion (about 1 sample radius in a 
time less than 1/24 sec). See Appendix C for n discussion of this. 


V. ANALYSIS OF MOTION OF SPHERE 


A. The Pressure Field 

The force on a sphere in a high -intensity , standing- wave acoustic field was 
first predicted by King 1 1] and has since been verified experimentally. (The exis- 
tence of radiation pressure in a sound field was originally predicted by Rayleigh, 
and the theory was developed in a series of 28 papers from 1902 to 1905.) The 
effect is one of second order in the fluid equations for the pressure in the sound 
wave. 


The basic equations are 


3 u . . 

Tf + ( H 


1, 


V ) u = -Vp 


»p 

3t 


+ V • (pu) = 0. 


where the viscosity term has been omitted from the first equation. This is accept- 
able for determining the pressure- velocity relationship but no* for predicting acoustic 
streaming (to be discussed later). 

For an irrotational fluid. King defines the velocity potential by 


u 


■Vi- 


and gets, in the simplest approximation. 


r , 2 , 1 ? 2 1- 

V * = ~2 — 

c at 


where c = /cip / dp is the speed of sound. Taking p to be the variation from ambient 

2 2 

pressure p Q , the result for p to second order in u /c is 


3* . 

P=!J oh + 


1 p o /a A 2 l 2 

2 T \3t j 2 f o- ’ 


where p is the ambient density. 
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King ignored any space dependence of <!■ other than that along the axis of 
sound source and reflector. If we generalize King's result to include a radial 
dependence, using cylindrical coordinates with cyclir. urical symmetry, we have 


<1> = A R(r) cos kz cos t, 


where A is constant. Since the standing waves in a St. Clair geneiator arc not 
confined to a well-defined cavity, we cannot treat the system by standard b -undary- 
value theory. We can, however, require that the function R(r) give a r idial depend- 
ence of the force on the sphere that corresponds to that which is observed. (We 
have done laboratory work on the force on a stationary sphere suspended from a 
microbalance, which can be combined with the work of Oran, et al. [2], to give a 
coherent picture.) 

Taking 


4’ = (r,z) cos a' t 


we get 

+ -> 1 2 r = 0. 

(Note that if there is radial dependence in .}> , k # w/c.) 

(k # .o/c) is apparent in data already available from Oran, et 
al. [2]. They measured the axial force in a sphere at room temperature and got 

force maxima at a separation of 1.41 cm, corresponding to a k value of 2.22 cm *, 

compared to the value of Wc (at 15 kHz) of 2.73 cm \ and giving y - 1.59 cm *. 

The radial equation given previously has a solution 


d*~R + 1 d R 

.2 r dr 
dr 


. 2 . u) 2 

where > -g - k . 
c 

Such an effect 


R(r) = J Q (>r) > 


where is the zeroth order Bessel function and > is to be determined experiment- 
ally. Using this result in the pressure equation, we get 


1 ' o 9 2 2 2 2 

p = -p u'A J (>r) cos kz sin u>t + =• -5 A J (>r) cos kz sin u >t 

OO MW O 

c 

1 2I22 2 22 2 T 2 

- g p A y Jj ( y r) cos kz + k J q (^,r) sin kz cos *ot , 
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With -g£ = - y.^ (yr>. We get tht "dc" pressure by averaging over one or more 
cycles : 

1 2 . 2 [" T 2 . ¥ 2 , .1 2, ^ 1 ,2.2, 2 . . _ 
p = ^ p© Y A J q (yr) - Jj (yr) cos kz + ^ p k A (yr> cos 2 kz. 

The last term, with J Q 1. is the result obtained by King, valid for yr << 1. 

If we take the force on the sphere to be given by 

F ' - (const) V Y, 

we can compare the radial and axial parts. This is equivalent to defining a 
potential -energy function V(r,z) with V(r,z) « p and taking 

V(r ,z) = V + 1 k + k r“ + .... 
o 2 z r 

We find k /k at the equilibrium point (near r = 0, z = rr/2k) to be given by 

Z I* 

k = 2K (2k 2 + Y 2 ) k 2 
z ' 

, 2, 2 
k r = K y k 

k z /k r = 2(2k 2 + y 2 )/y 2 . 

If we determine the periods of small radial oscillations, we have 

. ,, _ „ 2,_ 2 . . 1 T r 2 .2 „ . 2 

k z /k r a r /T z ’ Pving 2 — J - 1 Y = 2 k 

T z 

2 2 2 2 

and , of course, k + y = u» /c . Taking the data from SPAR VIII, we get 

T IT : 3. This gives 
r z fJ 

k = A y/2. 

With to/c = 1.09 cm * at 1575°C, we get 

k = 0.87 cm y = 0.66 cm 
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We can use these values to estimate the departure of the potential energy func- 
tion from Hooke’s law. There is obviously more to be done in the study of the 

5 

pressure function, including a reconciliation of King's result F = -j V dp/dz. 

Z 4 

Plots of the force and potential from the preceding calculations are given in Figures 
30 through 33. 


B. Orbit of the Sphere 

Consider the radial motion of the sphere. Since k r and k z are essentially 

constant, we can treat the axial motion and the motion in the radial plane as inde- 
pendent. We approximate V(r) by 

V(r) =|k p r 2 . 


With the substitutions 


u = 1/r 


du _ du fl 
u " dt - d0 0 


the equation for the orbit becomes [3] 


e - e = - 

0 


/ 


du 


/2mE 2mV 2 
u o • 7 - - 


2 . 


The substitution x = u gives an integral of the form 
x 

I = 


/ 


du 


/ — r 5 

/ a + bx + cx 


with the result 


1 

1 


2 2 


i/I* 


- 4B sin 2(6 - 0 


2(0 - 0 Q ) , 


where A = B = — 

r z 


mk 
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The energy equation is [1] 


E = l mr' 2 + - — f, + i k r 2 , 

2mr^ 

and so the turning points of the radial motion given by [4] 

,2 


E - - *1 + i k r 2 

* T 2 r 1,2 


or 


2mr 


£ - ^ + I /a 2 - 4B 

2 2 " 2 


1,2 


1,2 


The simplest form is obtained by choosing 0 q =• v/4, so that 


^ j /a 2 - 4B (cos 2 6 - sin 2 9) , 
r 

which is easily shown to be an ellipse with the center at the origin, major axis along 
the x axis, and semi-major and semi-minor axes and r ^ given by 


1 _ 1 A 1 /7 2~ 

~2 ” 2 A ‘ 7 ^ A 
r l 


4B 


"V 5 A + 1 /a 2 - «. 

r 2 


The period x of the motion follows from Kepler's second law [2]: 

2m, 


r - — 


where S is the area of the ellipse. 

„ /k 

_ 2tt _ / r 

Ul) = “7 / . 

r £ n 


We find 


The idealized motion of the sphere (no perturbations, no departure from 
Hooke's law force) is thus an elliptical orbit with angular velocity modified by 

axial oscillations with angular frequency a> = / k z /m . 

The problem is also separable in Cartesian coordinates, and we get for the 
time dependences of the motion , 

x = A^ cos (ul^t + 4>^) 


y = A 2 cos (uo r t + < t> 2 ) 
z = Aj cos (w 2 t + <t>g) . 


If the ellipse is oriented as given previously , the x and y equations are 


x = r^ cos at 


y = sin wt. 
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C. Perturbations of the Orbit 


1 . Time independ ent corr ecti ons to Hooke's law_. 

Bertrand proved [3] that the only central forces that give exactly closed orbits 
are the inverse -square force and the Hook's law force. Thus, any modification to 
the potential in Part B will give an orbit that is not closed. The time-independent 
modification takes the form of a fourth- power term in the expansion of the radial 
potential. For approximately circular orbits the effect of such a correction can be 

estimated. The correction to u h - is given by [3] 
u = u + a cos 30, 


where 3=2 for a Hooke's law force. One can then shdw [4] that 

3 2 = 3 + ^- ^- r I _ r , 

F dr r - r 
r o 

where r Q is the average value of r and F r is the radial force. 


From the expansion 


J 

o 


2 


(x) = 1 


1 2 
2 X 



4 

x - . . . 


. 12 2 
= 1 - 2 Y r 


16 


4 4 
Y r 


and the definition of a fourth-power coefficient, 


V(r) = V + k d r 2 
o r i 


2 bp 4 + 
4 




we see that b = ^ and 


dF 0 

r- ~5r - 1 ’ 2br2 + • • • 
r 


This gives 


6 2 = 4 - 2br o 2 = 4 - | Y 2 r o 2 


or 


V 6 1 = 2 - 3 ~ 0.04 
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for y = 0.66 cm 


and r 

o 


* r 

max 


The significance of this result is that the correction term causes the elliptical orbit 
to precess in the angular direction of the orbit. Qualitatively , this occurs because 
angular momentum is conserved; but the sphere stays near r mftx longer than it 

would for a pure Hooke's law force. The amount of precession, measured in the 
fraction of one complete cycle of precession for one cycle of the elliptical motion, is 
V 8. That is, after one cycle of 0 in a space-fixed reference frame, the sphere has 
not quite returned to its original radial position. Our estimate shows that the pre- 
cession could be up to '200° in 14 cycles, or slightly more them ^ of a revolution. 

During the total motion of the sphere, the amplitude could thus go from a minimum 
to a maximum and back, or vice versa. 

2. Time-dependent force constant . 

The data on the position of the sphere indicate a time variation of the force 

constants k and k in such a way that k / k ~ constant. This indicates a decrease 
z r J z r 

in the dc pressure in the sound field, the explanation for which is discussed else- 
where. We consider here the effect of this time -dependence of k and k on the 
motion . z r 


The equations of motion retain their same form and are modified only by the 
fact that the force constants are time dependent: 


mf 



yt) r 


mz = - k (t)z 
z 


with l, the angular momentum in the plane of the orbit, still constant. 

Starting at t = 5 sec = t , we find from the data that k (t) makes an abrupt 

o z 

change of about 50 percent during the first 5 sec and a gradual change of about 
40 percent during the next 60 sec, after which the motion becomes quite erratic. 
Let us consider the middle portion during which the change is relatively slow. 

If k z (t) changes sufficiently slowly, we can represent the motion by 


A(t) = Aj(t) cos [u(t) t + <J>(t)’ 

(adiabatic approximation). This approximation will be good if Aw/u)<<1 for one 
cycle, a condition which holds reasonably well during the period from t = 10 sec to 
t = 70 sec. Einstein proved [3] that in the adiabatic approximation the ratio of the 
mean energy to the frequency remains constant. The problem was earlier solved by 
Rayleigh as the starting point of his study of radiation pressure (1902). 



Einstein's result also follows from the Boltzmann -Ehrenfest theory of adiabatic in- 
variance (1913). The energy in this approximation is (considering the axial motion 
only) 

p - 1 m 2 a 2 

E z “ 2 m u A 3 ’ 

2 1 

so we get loAj = const or Aj “ - . The change in w is from ~ 3.5/sec to 

* / cr 

o >2 = 2. 6 /sec. This gives about a 20 percent increase in A^. The actual behavior 
of Ag is quite different from this, indicating that some other perturbations are 
dominating the amplitude but not significantly affecting the period. 

3. Numerical analysis. 

The preceding analysis indicates that the potential function describing sample 
motion in the acoustic levitation system is time dependent and requires terms in 
the radial equation higher than second order. For the radial potential function we 
use the form: 


V(r,t) = K q (^-rit) Ar 4 ) • 

From the flight data K is determined to be approximately 1.5 in cgs units. B is 

0 - 3 - 1 

estimated from the time behavior of the period to be approximately 8.6 x 10 sec . 

A is obtained by applying the solution of the time -independent correction (discussed 

in Section V.C.l) to the geometry of the flight acoustic levitator. A is estimated to 

. 2 

have the value of -0.33 cm . Other parameters are the sample mass (0.6 gm) and 

an estimate of initial orbital angular momentum of the sample (<0.2 gm cm sec *). 
The estimate of angular momentum comes from assumption of near-circular orbit and 
measurements of sample mass and orbital period. 

The force derived from the preceding radial potential function is: 

1 2 3 

F(r,t) = -K o (j-pgt) (2r + 4 Ar J ) . 


The equation 


2 2 

m ; - L_ = - Kq ( ir J- R ) (2r = rAr 3 ) 
mr 

was solved numerically, using the parameters shown previously, to yield a sample 
orbit. Figure 34 shows the first 2i cycles of the computed orbit. The important 
qualitative feature that can be seen here is that an orbit approximating a precessing 
ellipse has been obtained. The integration has been carried out for 90 sec, and a 
projection of the resulting orbit onto a line in the orbital plane is shown in Figure 35. 
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This projection of the orbit reveals some features seen in the flight data. First, 
there is a secular amplitude variation reminiscent of the variation seen in the flight 
data (Figure 26). This is the result of observing the orbit from the orbital plane 
and seeing projections of the orbit which vary as the orbit precesses. We also see 
an increase in the orbital period with time, as is noted in the flight data (see Figure 
27). Figure 36 is another projection of the radial component of motion with the same 

_2 

input parameters as previously, except that A = -0.5 cm . The result is qualita- 
tively the same as before, but the increase in the orbital period is faster than 
before. The change in A represents a shallower potential well, as indicated in 
Figure 2. 

2 

If A = 0 (that is, the potential is restricted to an r dependence in its spatial 
variation), the first few cycles of the resulting sample orbit will be as shown by 
Figure 37. We see that the precession of the orbit is now very slight. The fairly 
rapid precession implied by the radial amplitude variations then dictates the need 
for spatid variations of higher than second order. 

4. Drag. 

There will be drag forces on the sphere because of its motion relative to the 
air, which itself also has a time-independent flow known as acoustic streaming. The 
force is usually represented by the equation 

F D = \ p o A y2c D 

= (3.02 x 10 5 ) v 2 Cp (cgs units), 


where Fjj is the drag force, A is the cross-sectional area of the sphere, v is the 

speed of the sphere relative to the air, and C D is the drag coefficient [5]. The 

drag coef'-cient depends only on the Reynolds number R (for uniform flow without 
boundaries) , defined by 


R 


Dv 

v 


where D is the diameter of the sphere, v is the kinematic viscosity, and v = y /p 
where y is the (dynamic) viscosity. For our sphere at 1575° 


R = 0.204 v 


for v in cm /sec. C D has been determined for a broad range of values of R and is 
plotted, e.g. , in Reference 5. 

We first estimate the drag force which would exist in still air due to motion 
of the sphere. A typical maximum speed is obtained from the amplitude and 
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Figure 36. Projection of calculated orbit onto a line. 
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frequency data. Taking radial and axial amplitudes of 0.5 cm and 0.25 cm and angular 
frequencies of 1.5/sec and 5/sec, respectively, as near maximum values, we get a 
maximum speed of 1.5 cm /sec. (We will see later that acoustic streaming numbers are 

_ - 3 

larger.) This gives R ~ 0.3, C n 100 and F - 5 x 10 dynes. Typical value.' 

u o 2 

during the motion will be much smaller. (Stokes' law gives F q = 3.6 x 10 v.) 

Thus it appears that we can ignore the motion of the sphere itself and consider 
only acoustic streaming when estimating drag effects. 

5. Spin . 

Since the sphere was determined to be spinning during part , if not most of the 
motion , we should also consider the force on it due to its spin (Magnus effect ; see 
Reference 5, p. 71). The force per unit length T^/SL on a rotating cylinder moving 

at speed v is given by [ 5] 


f n = p r v. 


where Y is the circulation of the cylinder [5]. We can calculate the equivalent effect on 
a sphere as 


8 „ 3 

Vi' p o vu "3 

= 5. 12 x 10 ^ Vf, (egs units) , 
s 

where a is the radius of the sphere and is spin angular velocity. For motion due to 

s 

the sphere alone, using v = 1.5 cm /sec as before and an estimated value for w based 

s 

on a (guessed) rotational frequency of 10 Hz, we get 
F. Z 3 x 10 ^ dynes. 


Thus we could increase io by two orders of magnitude and still ignore F . We will 
reconsider spin, however, in connection with acoustic streaming. 

6. Acoustic streaming . 

High -intensity sound waves in air are accompanied by time- independent (dc) 
flows known as acoustic streaming (see, e.g., Reference 6 and 7). These flows occur 
in a free nonuniform sound field, near objects in the field, and near oscillating bodies. 
They are intimately related to the subject of radiation pressure which gives the force 
used for positioning of the sphere. They were first noted by Faraday (1831) and were 
studied by Rayleigh in connection with radiation pressure [81. 


For streaming around an object in a sound field, there is a dc boundary layer 
which separates two flows. The inner flow is know as Schlichting streaming and the 
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outer one hs Rayleigh streaming. Flow in e free non uniform sound field is known as 
Eckart streaming. For all three *vpes, stability is governed by viscous forces. 

The nature of streaming 1 a standing wave field is known qualitatively and 
quantitatively at low intensities, and some work has been done at high intensities. 

An estimate of the magnitude of the dc velocity can be obtained from the measure- 
ments of Borisov and Statnikov [9] at high intensities. They used three different 
methods to measure the axial component of the dc flow velocity in a circular tube and 
reported the maximum values versus dc pressure, giving 

- , . -2 

u - b p , 

- 3 2 2 

with b = 2.2 x 10 (cm /sec) /(dyne /cm ) . The curve was determined for sound 
pressure levels between 163 and 168 db, with measurements at lower SPL giving 
values smaller than those predicted by the curve. 

The basic relationships are: 
sound pressure level: 

SPL = 20 log 1A (P / p ), 

°10 rms *s 


where P is the root -mean -square pressure in the sound wave and p is the stand 
rm 8 _ 42 ® 

ard reference pressure d = 2 x 10 1 dyne /cm . 

s 


dc pressure: 

f : p 


2 *> 

/ 2p c , 
rms o 


2 

where p q c = y P q for ideal gas, with > = ratio of heat capacities of air and P q = 
mean reference pressure = mean pressure in payload. Thus, 


u « 


P 

rms 


4 


3 

according to Borisov and Statnikov. For a sample density p^ of 4.43 gm/cm 

calculate the approximate value of P using an equation given by King [1], 
result gives rms 


we can 
King's 


P 

rms 




5 p 


^ R n 


» 


where n is the angular frequency of axial oscillations of a sphere in the sound field. 

We have checked the relation of u to p given by Borisov and Statnikov in the 
laboratory at room temperature by photographing spheres of uniform density (50 pm) 
in the sound field at measured intensity . using a well-defined laser beam for 
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illumination. The results agree reasonably well. For a determination of u at 1575°C, 
however, we must also take account of the variation of kinematic viscosity with tem- 
perature. Reynolds' principle of similarity (see, e.g., Reference 5) gives u a v. 

For the same value of P , then, the expected value of u will be greater by the 

rms 

ratio of kinematic viscosities: 

v'/v = 3.12/0.15 = 21. 


The results , using the formulae of King and of Borisov and Statnikov with the 
preceding factor, give the values shown in Table 1. The values of n are the co 
values quoted in Section V.C.2. 

TABLE 1. ESTIMATED MAXIMUM VALUES OF ACOUSTIC STREAMING VELOCITY 


n(sec *) 

2 

P rms (dyne/cm J 

SPL(db) 

2 

p( dyne /cm ) 

u(cm/sec) 

3.5 

6.4 x 10 3 

150 

14.7 

10 

2.6 

4.8 x 10 3 

148 

8.1 

3 


From the discussion in Section V.C.3, we find that for u = 10 cm /sec, R = 2; 
this is still within the range where Stokes' law holds (force a velocity) , and we find 

- 2 

that typical forces on the sphere due to acoustic streaming are less than - 4 x 10 
dynes . 

Although our data on acoustic streaming are not yet definitive , it appears that 
the highest streaming velocities occur at the greatest distances from the equilibrium 
position . within the range of motion of the sphere. This is, of course, also where 
the restoring force is greatest. By comparison, the restoung force is - 1 dyne. 
Thus, while streaming forces are not of the same order of magnitude, they are not 
negligible in determining the shape of the position versus time curves and appear to 
be of about the same relative importance in this respect (shape, rather than period) 
as the precession and the time variation of the force constants discussed previously. 

It should be pointed out here that one cannot determine the streaming velocities 
from the equations in Section V.A. , since they are produced by the stress term [7J, 

2 

y V u, which was omitted from these equations in a calculation of the pressure. 

Whether we could make accurate predictions of the shape of the position versus 
time curves by detailed studies of the streaming is relatively unlikely. However, we 
believe that a knowledge of the direction and magnitude of the streaming velocities in 
high -intensity sound fields will be important for implementation of acoustic levitation 
devices. Certainly streaming is primarily responsible for the shearing forces which 
cause samples to spin, and we know that the spin can be affected by adjusting the 
angle and/or the position of the reflector. 
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VI. SUMMARY AND CONCLUSIONS 


A. Summary 

In this report we have given a description of the single -axis acoustic levitation 
experiment on SPAR-VIII. We reviewed the primary data, mainly a film of the 
sample motion, and various graphical and digital representations of these data. 

Having presented the data, we then identified the phenomena which, in our 
judgment, could or should be explained. This led to an analysis of the forces acting 
on the system, an estimate of their relative magnitudes based on a theoretical model, 
and a study of the effects of those forces which seem to be significant. 


B . Conclusions 

Here we list, briefly, the conclusions which we have drawn based on our anal- 
ysis, as discussed in more detail in the body of the report. 

1) There are no foi*ces due to anything external to the experiment package 
itself which had any significant effect on the motion of the sample. 

2) The shift in the axial equilibrium position was caused by the cooling of the 
air on injection of the sample and its subsequent return to operating tem- 
perature as the experiment progressed (see Appendix D) . 

3) The shift in the radial equilibrium position was due to a combination of an 
asymmetric distortion of the acoustic field caused by the injection device 
and the temperature change and axial shift discussed above (see Appendix 
D). 

4) The gradual increase in the radial and axial periods (with constant ratio) 
was due to a lowering of the sound pressure level. This resulted from 
increased water vapor due to evaporation from the injector and sample as 
they heated up. Conclusions 2, 3, and 4 are drawn pending further 
experimental confirmation from Inter sonics, Inc. 

5) The gross features of the motion were what one would expect from the 
theoretical model which was presented. They consisted of an elliptical 
motion in the radial plane coupled with axial oscillations. 

6) The modifications in the amplitude are due to a combination of effects. In 
the absence of three-dimensional data, these effects cannot be isolated. 

The primary ones are coupling between the rsdial and axial oscillations due 
to the nature of the pressure field, elliptical precession due to the fact 
that the radial potential is not purely simple harmonic, and drag due 
primarily to acoustic streaming. 

7) The apparently random changes in shape of the position vers us -time plots 
are due primarily to acoustic streaming. We draw this conclusion based on 
our estimated magnitude of the streaming velocity and the elimination of 
other effects from consideration. 
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8) The finul loss of the sample was due to a combination of weakening of the 
sound field and i.ie attractive force of acoustic streaming when the sample 
neared the cage. 

9) The final, rapid motion of the sample was a wetting effect after the i olten 
sample contacted the cage. 


C . Recommendations 

The following recommendations are made: 

1) Reduce the input of water vapor to the sound field during and after injec- 
tion. This could be accomplished by purging the package with dry air 
prior to flight and sealing it. or by preheating the sample and injector 
during flight to 800 degrees C. 

2) Expand the cage so that it is outside the range of the restoring force of 
the sound field. This will insure that the sample will not be drawn to the 
cage by acoustic streaming unless it is first lost from the sound field. 

3) Provide a second view of the sample so that three- dimensional motion can 
be reconstructed. This should be accomplished by a mirror arrangement 
and composite photographs using a single camera. 

4) Use color film to increase contrast. 

5) Check the injector design in the laboratory to assure that it does not 
significantly distort the sound field, or plan for the distortion. 

6) Test the experiment package on one or more KC-135 flights exactly as it is 
to be flown on a SPAR flight. This test should include an analysis of the 
motion similar to the one in this report, so ns to give an accurate predic- 
tion of the gross features over the longer period of time on a SPAR flight . 
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APPENDIX A 


The sample used in the flight experiment was reported by Dr. Ralph Happe of 
Rockwell to be composed of gallia (41 percent), calcia (36 percent, and silica (23 
percent), where percentages are mole-percentages. Mass was reported to be 
0.595 gm. Measurements of sample diameter and density were unavailable. These 
were computed as shown below from densities of the components as reported in the 
Handbook of Chemistry and Physics, 

densities: 

3 

gallia- 5. 88 gm/cm 

3 

calcia-3.30 gm/cm 

3 

silica- 2.20 gm/cm 


Define: N , N , N as the mole fractions of gallia, calcia, and silica, respec- 

A d L 

tively, and M^, Mg, and N^. as the molecular weights. Similarly, define C^, Cg, 
and Cg to be the weight fractions of the respective components. Then: 

N a M. 

r _ A A 

A N . M A + N 


A A 


B M B + N C 


c. = 


(0 . 41) ( 188) 


_ 77.08 _ 


A ( 0 . 41) ( 188) + (0.36) (56) + (0.23)(60) 111.04 


0.694 


. (Q .36)(56) _ 

B 111.04 C - 182 


o _ (0.23)(60) _ i0 . 

C C " 111.04 ~ ° 121 * 


The density cf the sample is given by 


?A + S + ^C 


P ABC P A °B 


where p. , p , and o are the densities of gallia, calcia, and silica, respectively. 
Ad l' 


°ABC 


.694 , .182 .124 

5.88 ‘ 3.30 2.29 


\bc = 4 - 43 e m/cm ' 


From the sample density and weight, the sample radius is calculated to be 
0.318 cm, or 0.250 in. 
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APPENDIX B 


The projection of sample trajectory onto the y-z plane is shown in the series 
of plots in Figures B-l through B-20. Each plot represents position data from 100 
consecutive motion picture frames, or approximately 4.17 sec of data. Each plot, 
then, contains approximately two-thirds of a radial period and approximately two 
axial periods (on the average) . 



AXIAL (INCHES) 



AXIAL (INCHES) 

B-2. Sample trajectory for time t = 8.38 to 12.50 sec. 
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AXIAL OnotMM) 

B-ll. Sample trajectory for time t = 42.39 to 46.52 sec. 



AXIAL IMm| 

B-12. Sample trajectory for time t = 46.56 to 50.69 sec. 
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APPENDIX C 


i 


The sample is observed to move a distance approximately equal to its radius in 
a time less than 1/24 sec; that is, the sample moved that distance between successive 
motion picture frames with the camera operating at 24 frames per sec. After this 
motion, the sample position remained fixed. The interpretation placed on this is that 
the sample contacted the wire cage and, because of capillary forces, flowed around a 
wire and remained centered on the wire because of these forces. 

The force required to uniformly accelerate a mass over a distance d in a time t 
is given by : 



for this sample: 


f 


(0. 595 gm. ) 


(2)(0. 318 cm) 
U/24 SeC > 2 

k 


f = 218 dynes. 


The acceleration of the sample is: 


f = 218 djnes = on:/se< ,2. 

m 0.595 gm 


In terms of g, this is: 


366 n 

a = 980 * " 0 37 «• 


This acceleration is too high to come from any available source except surface tension. 
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APPENDIX D 


An approximation of the drift in the average radial position of the sample as a 
function of time is given by: 


R = 475 + 285 e 


-0.035 


9 


where R is the sample average position in thousandths of an inch. A plot of this 
function is shown in figure D-l. Similarly, an approximation for the average axial 
position as a function of time is given by: 


A 


800 - 


350 e' 0,0553 * 


This is shown in Figure D-2. 

If we make linear approximations to the radial and axial force constants (from 
Figures 24 and 29) , we may write 


K d = 1.10 - 0.0095 t 

K 


and 


K. = 11.0 - 0.095 t, 

A 

2 

where K R and have units of dynes /cm or gm/sec . 

If we assume that the drift in the sample average position is due to some 
external force superimposed on the acoustic field, we may use the expressions for 
R, A, K„ , and K. to compute these forces as functions of time. They may be 

A A 

expressed as : 


Fr = k r ar 

and 


F A = k a aa. 
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o«ition (inoh##) 






Converting R and A to cgs units, we have: 


F r = [2.54 x 10~ 3 ][1.1 - 0.0095 t ] 1 285 (1 - e' 0,035t )] 


and 


F a = [2.54 x 10' 3 ][11.0 - 0.0 95t ] [350 ( e "°- 0553t -1)] . 

These forces are plotted in Figures D-3 and D-4, respectively. 

To compare the magnitudes of these forces to others operating in the system, 
it is convenient to normalize them with respect to gravitational forces. In terms of 
g, the accelerations a mass m would undergo because of these forces are: 



JU L 


980 cm /sec 


and 


. A .i* 


Jl g_ 


l 980 cm /sec' 


For the sample with m = 0.595 gm, we have: 


a R = (0.0017 F r ) g 


and 


a A = (0.0017 F a ) g . 

-4 -3 

From Figures D-3 and D-4, these have values of approximately 10 g and 10 g, respec- 
tively. Since these values are so high (1 to 2 orders of magnitude greater than 
measured values), we can conclude that the observed drifts in average sample posi- 
tions do not arise from external forces. They must be due to changes in th«* 
acoustic field itself. 
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Figure D-4. Apparent axial force. 
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